clear all; clc; close all;
N = 10000;
sigma = 2;
x = 10+(22-10)*rand(1, N); %1пункт
%x = sigma*randn(1, N)+16; %2пункт
x1 = sigma*randn(1, N)+16*0;
x2 = sigma*randn(1, N)+16*0;% 0-Релей, 1 - РАйс
%x = x1.^2 + x2.^2; %3 пункт
x = sqrt(x1.^2 + x2.^2);




hist(x, 30)
break

plot(xk,  p);
grid on;

xlabel('x'); ylabel('p(x)');

 